Energy landscape and phase transitions in the self-gravitating ring model 
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We apply a recently proposed criterion for the existence of phase transitions, which is based 
. on the properties of the saddles of the energy landscape, to a simplified model of a system with 

gravitational interactions, referred to as the self-gravitating ring model. We show analytically that 
the criterion correctly singles out the phase transition between a homogeneous and a clustered phase 
and also suggests the presence of another phase transition, not previously known. On the basis of 
^ ■ the properties of the energy landscape we conjecture on the nature of the latter transition. 
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Phase transitions (PTs) are abrupt changes undergone by the collective properties of a system at a given value 
£j i of a continuously varying parameter, like temperature or pressure. Common examples are the melting of ice or 
the spontaneous magnetization of a ferromagnet. From a theoretical point of view, PTs are usually defined as 
£h nonanalyticities in thermodynamic functions [25j. 

Our current understanding of many aspects of PTs is remarkable This notwithstanding, there still are some 
open problems that should not be overlooked. One of such problems is that we do not know which features of the 
' Hamiltonian of a system induce a PT. The interaction potential energy V(q±, . . . , <?at) among the microscopic degrees 
. of freedom contains, in principle, all the information on the behavior of a system, at least when the Hamiltonian 
is separable and the kinetic energy term depends only on the momenta, that is for most classical many-particle 
systems (iij. Given V, however, we typically do not know what to look for in order to say whether the system has 
or not a PT. Some necessary conditions on V are known, like those expressed by the Mermin- Wagner theorem Q 
excluding PTs for certain classes of potentials with a continuous symmetry or by the Van Hove theorem Q (later 
generalized by Cuesta and Sanchez [4j) excluding PTs for suitable one-dimensional systems, but almost nothing is 
O . known about sufficient conditions V must fulfill for a PT to take place. 

A natural framework for approaching this problem is energy landscape theory [5[ . The energy landscape (or more 
H ■ precisely Potential Energy Landscapes, PEL) of a classical system is the graph of the potential energy function 
£> ' V : Tjv i—* K, where Tn is the configuration space of coordinates q%, . . . ,q^, and one searches for those features of 
the graph which may be associated with the relevant physical properties of the system. This has been fruitful in 
the theory of glasses and disordered systems, whose complex behaviour may be traced back to the presence of many 
near-equivalent minima in the PEL [fjj , as well as in protein folding theory, where the ability of a protein to fold in a 
unique native configuration is ascribed to the overall funnel shape of the PEL 0] , with the native state at its bottom. 
Not only minima, but also other stationary points of the PEL, i.e., configurations such that dV = 0, play a key role 
' in energy landscape theory and are often referred to as saddles. Minima (saddles of order zero), or better their basins 
of attraction, identify the regions ("valleys") of r^v which are visited with a high probability, while saddles of higher 
order are associated to transition states between different valleys in the PEL. Saddles of higher energy may become 
• i-H , accessible upon increasing the energy or the temperature, so that different valleys merge and the shape (topology) 
of the effectively accessible configuration space changes. When such a change is sufficiently st rong , it may entail a 
PT; this is the basic idea of the topological hypothesis about PTs, first suggested in Q (see [^, EI EH for reviews). 
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Indeed, Morse theory [12| can be used to show that there is a one-to-one correspondence between saddles of the PEL 
and topology changes in the manifolds M v ={g£ Fjv|F(g) < Nv}. In addition to many calculations carried on for 
model systems (see refs. in 13, two results have given a solid basis to the topological approach to PTs. First, 



Franzosi and Pettini [3] proved that under rather general conditions topology changes in the M v 's (thus stationary 
points in the PEL) are a necessary condition for a PT to occur. Second, Kastner, Schnetz and Schreiber (KSS) [14[ 
proved that the microcanonical entropy of a system with a finite number N of degrees of freedom is nonanalytic at 
every stationary value of the potential, so that there is a one-to-one relation between topology changes of the M v 's, 
stationary points in the PEL and microcanonical singularities. However, there are many stationary values of V (their 
number being typically proportional to e N ); moreover, they usually become dense in all the interval of values of V 
as N — > oo and, in the same limit, the majority (if not all) of the microcanonical nonanalyticities become irrelevant 
since they show up as singularities of N-th order derivatives of the entropy [3, EH . This suggests that the definition 
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of a PT as a nonanalyticity of a thermodynamic function is unsatisfactory in the microcanonical ensemble and should 
be restricted to singularities surviving to the thermodynamic limit; at the same time, most of the saddles of the PEL 
(and of the corresponding topology changes of the M„'s) turn out to be unrelated to PTs. 

A major step in understanding which are, among all the saddles in the PEL, those that may induce a nonanalyticity 
in the entropy which survives in the limit N — > oo and that may then be related to thermodynamic PTs was provided 



again by KSS [14 . [16| . They put forward a discriminating criterion whose content is essentially that such saddles must 
become sufficiently "flat" as N — ► oo. The KSS criterion, as well as most (if not all) of the other tools proposed within 
the topological approach to PTs, has been successfully applied 14, 1(| until now only to models that are (i) exactly 



solvable and (ii) whose energy landscape is completely characterized, i.e., essentially all the stationary points and 
stationary values can be analytically calculated. Such models, although clearly important as a conceptual benchmark, 
are highly nongeneric. 

The aim of the present Letter is to analyze further the potentialities of the KSS criterion by applying it to a model 
system that, contrary to those previously studied, is not exactly solvable and whose energy landscape is not completely 
characterized analytically. Nonetheless, the criterion is able to single out the stationary points associated to a PT; it 
correctly detects the homogeneous-clustered PT and also suggests the possibility of another PT, previoulsy unknown. 
Moreover, in spite of the fact that the model is not exactly solvable so that all the information on its equilibrium 
behavior must be obtained numerically, our analysis has been carried out analytically. Our results then show that 
the KSS criterion can be effectively applied to generic systems, more complex than those previously addressed in the 
topological approach to PTs, and that it may even provide new insights and exact results otherwise unattainable. 

The KSS criterion. Let us briefly discuss the KSS criterion (for details see Consider a classical Hamiltonian 

of the form 

1 N 

w(p,«) = 3!^ + % «f) w 
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and assume that the stationary points of V are isolated and that their overall number grows at most exponentially 27 
with TV. Then the microcanonical entropy per degree of freedom s(e) can be nonanalytic in the limit N — * 00 at e = e c 
only if the following two conditions are met. First, there is a sequence {q^} N _ 1 of stationary points of V whose 
corresponding stationary values converge to v c = (v)(e c ), where v = V/N and (v) denotes its statistical expectation 
value, i.e., 

Jim v (q?) = v c . (2) 
Second, the Hessian matrix Hy of the potential V computed on the stationary configurations is such that 
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lim |det[H v (q?)]\ =° • ( 3 ) 



Since the eigenvalues of Hy can be seen as curvatures of the PEL, Eq. Q means that the saddles become asymptotically 
"flat". We note that to check whether Eqs. @ and arc satisfied we do not need to know all the saddles of the 
PEL; it is sufficient to determine the "right" ones. This is a big difference with respect to other criteria previously 
proposed, which made use of topological invariants like the Euler characteristic whose calculation requires to know 
all the saddles of the PEL (and their order) @,[llj]. 

The model. The model we studied is the Self-Gravitating Ring (SGR), first introduced in as a simplified model 
of a self-gravitating system, ft is a model of N points of unitary mass moving on a circle of unitary radius and 
mutually interacting via gravitational forces, regularized at short distances. Its Hamiltonian is of the form |T|) with 
potential 
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V = — V (4) 

2NV2 ^ v/l-cosfai-fcO + a ' 

where qi € (— n, 7r], i — 1, . . . , N, are the angles giving the position of the «-th particle on the ring and a > is the 
softening parameter regularizing the potential for (qi — qj) — > 0. The factor in Eq. ([4]) ensures extensivity according 
to the Kac prescription, although the model is non additive, the interactions being long-range (l9| . 



Numerical results (iV-body simulations 18] as well as variational calculations in the N — > 00 limit [20() for the 
microcanonical thermodynamics show that the system has a PT separating a homogeneous high-energy phase from a 
clustered phase where the rotational symmetry of the potential is spontaneously broken. If the softening parameter a 
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is sufficiently small the clustered phase contains an energy interval with negative specific heat, there is nonequivalence 
between the microcanonical and the canonical ensemble and the qualitative behavior of the system is very similar 
to that of three-dimensional self-gravitating models. In the opposite limit a — ► oo SGR becomes equivalent to the 
ferromagnetic HMF model [2l[ . Assuming that in the high-energy phase the particles are homogeneously distributed 
on the ring, as suggested by the numerical results, one can calculate the expectation value of the potential energy 
density at the transition in the limit N — > oo, obtaining 
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where K(x) = f* /a ffl/y/l — x sin $ is the complete elliptic integral of the first kind. 

Saddles of the landscape and PTs in the SGR model. Let us now study the PEL of the SGR. First of all let us 
note that the stationary points of the potential ((J) are not isolated, due to its rotational invariance. However, this 
difficulty can be circumvented by fixing the value of one of the g's, which has an irrelevant effect on the thermodynamic 
functions as N — > oo and only fixes the position of the center of the cluster in the broken-symmetry phase: from now 
on we assume qi = 0. The stationary points of V are the solutions of the form (q\ = 0, q 2 , . . . , q^) of the N coupled 
nonlinear equations VV^ = 0, i.e., 
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with k — 1, . . . , AT. Physically the above equations mean that the force acting on each particle is radial. There are 
at least two classes of solutions that can be easily found. First of all there are the solutions we shall refer to as 0-ir 
saddles, where N v particles are in q = it and the others are in q — 0, with < N v < N — 1. We shall denote such 
configurations as q n7r , where iij = N^/N. The other class of solutions of Eqs. © that is easily found is that of 
configurations where the particles are in the p vertices of a regular polygon, with the same number r of particles in 
each vertex; from symmetry considerations one sees that the force can be only radial. We shall refer to the latter 
solutions as polygonal saddles and we shall denote them as q p<r , with N = pr. There are also many other saddles of 
the SGR potential that do not belong to any of these two classes. For instance, consider three angles 0, 7 and 5, and 
put Nq particles in 0, iV 7 in 7 and Ng in 5. It can be shown that for almost any value of 7 and 5 such that < 7 < 7r, 
— 7r < 5 < and < 6 — 7 < ir one can choose sufficiently large N, No, 7V 7 and Ng such that this configuration is 
stationary; moreover, one can find arguments suggesting that also more complex stationary configurations exist [22j ]. 
The polygonal and 0-ir saddles, then, arc not the only saddles of the PEL of the SGR model. 

This does not prevent the KSS criterion to be effectively applied to this model. Let us calculate the stationary 
values v — V/N corresponding to the 0-tt and polygonal saddles, respectively. For q niv we have 



v(n^) = - 
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The maximum of v(nn) is attained for = 1/2 and the minimum, corresponding to n n = 0, is also the absolute 
minimum of the potential. Hence 
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< v(n„) < - 
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and the values of w(rt T ) become dense in the above interval as N — * 00. For q p ^ r we have that the stationary values 
depend only on the number of vertices p: 
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The function v(p) is monotonously increasing with p and 



v(p) > 
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so that v(p) > v(n n ) \fp, n w . Moreover, a simple calculation shows that limp—^ v(p) — v c , where v c is given by 
Eq. ([5]) and is also the upper bound of the potential energy per particle. We have thus two nontrivial results: (i) 
since as N — > oo the distance between two successive stationary values tends to zero in both cases of v(n n ) and v(p), 
although 0-7r and polygonal saddles are not the only stationary values of V their stationary values do encompass all 
the available values of the potential energy of the model; (m) the sequence of stationary values v(p) converges to the 
critical potential energy v c of the PT between the homogeneous and the clustered phase as p — + oo. 

The last result suggests to investigate whether the KSS criterion is satisfied for v = v c , i.e., whether a sequence of 
saddles satisfying Eqs. ([2]) and ((3]) exists, with v c given by Eq. ([5|). To this end, consider the sequence of polygonal 
saddles with one particle in each vertex, qN,i, where N is prime; Eq. ((2]) is satisfied — see item (ii) above. It remains 
to show that also Eq. ((3|) is satisfied. Define 
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then the diagonal elements of the Hessian matrix are [Hy (lN,x)]uu — ~ X^o 1 f k ano - the off-diagonal ones are 
[Hy {<lN,i)]if. = /fe-Zj so that the Hessian calculated in qN,x is a circulant matrix. Using the Hadamard inequality [23} 
to obtain an upper bound to the absolute value of the determinant of a matrix as the product of the Euclidean norms 
of its rows and observing that in a circulant matrix all the rows have the same norm, after some algebra we can write 



lim |detpv (qN,i)]\ 1/N < lim 
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The r.h.s. of the above equation can be written as an integral, so that Eq. ([3|) is satisfied if 
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the above result can be proved by explicit integration, so that the KSS criterion is satisfied for v = v c . 

The properties of the circulant matrices [24| can also be used to prove that Eq. (J3j) is not satisfied for all the 
stationary values v(p) ^ v c (details may be found in 22]). Hence, as far as the stationary values of V associated to 
polygonal saddles are concerned, we have shown analytically that the KSS criterion is satisfied if and only if v = v c . 
Since the system is not exactly solvable, this is, to the best of our knowledge, the only analytical indication of the 
presence of a PT between a homogeneous and a clustered phase in the SGR model. 

To complete the analysis of the whole range of potential energy values we now turn to the 0-n saddles. Since 
the stationary values ([7]) depend only on = N^/N, let us consider sequences of saddles q n7r with varying N at 
fixed n-x- The Hessian calculated on such saddles can be written as Hy (q n „) = O + A, where A has rank 2 and 
D = diag(di, ...,d%) with 
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here Ai = 2V2a 3 and A2 
contribute to |det [Hy (q n7T 



2y/2{2 



Using such decomposition one can prove that only the elements of 
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where 0(71^) = n^/Ai — (1 — n n )/\2 and b(n n ) = (1— n w )/Xi — n 7r /\2. The quantity in Eq. (|15[) vanishes if and only 
if n-n — n% or n T = 1 — where n% — a 3 / 2 / [(2 + a) 3 / 2 + a 3 / 2 ] and v{n%) = v(l — n%) = v' c , with 



4 + a[6 + a(5 + 2a)] 
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The KSS criterion is then satisfied not only at v c given by Eq. ([5]) but also at v' c given by Eq. (fT6|) . Is there a PT 
in the SGR model at e' c such that (v)(e' c ) = v' c , as suggest ed by the KSS criterion? We do not have a final answer 
yet. Previous numerical studies of the SGR model [18|, l20( have not detected such a PT. However, we note that v' c 
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is extremely close to the absolute minimum of the potential for small values of a, where the SGR model is close 
to a "true" self-gravitating system. Conversely, v' c ~ v c for large values of a, where it behaves like a mean-field 
ferromagnet. Only for a ~ 1 one has that v' n i s clearly separated from both the mimimum and v c , but these values of 
a have not been thoroughly studied in [n| H3] • Moreover, the analysis of a related model [2§| suggest that such a PT 
may indeed exist but its effect on the thermodynamic quantities may be weak, so that it may not be easy to detect 
numerically. If such a PT exists, which phases does it separate? An answer may be suggested again by an analysis of 
the PEL. Computing the number of negative eigenvalues of the Hessian one can show that saddles q n ^ with n n < n£ 
are proper saddles, i.e., saddles of order greater than zero, while saddles with n v > n% are minima. Minima are visited 
with high probability, at variance with higher-order saddles. Hence we conjecture that for (v) < v' c the equilibrium 
phase is such that the fraction of particles that may cross q — n, i.e., visit the whole circle, is zero, while it becomes 
nonzero as (v) > v' c , i.e., at e > e' c . 

Summary. We have applied a recently proposed criterion for PTs, referred to as the KSS criterion, to the SGR 
model which, at variance with the few models the criterion had previously been applied to, is not exactly solvable. 
We have shown analytically that the criterion correctly singles out the PT between a homogeneous and a clustered 
phase. We have also shown that the criterion indicates the possible presence of another PT, not previously known, 
and we have conjectured on its nature. Our results suggest that the criterion may be effectively applied to even more 
complex systems, if one can guess which stationary configurations can be related to PTs. 
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